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Abstract 

In the first part, the induced vacuum spin around an Aharonov-Bohm flux string in 
massless three-dimensional QED is computed explicitly and the result is shown to agree 
with a general index theorem. A previous observation in the literature, that the presence 
of induced vacuum quantum numbers which are not periodic in the flux make an integral- 
flux AB string visible, is reinforced. In the second part, a recent discussion of chiral 
symmetry breaking by external magnetic fields in parity invariant QED3 and its relation 
to the induced spin in parity non-invariant QED 3 is further elaborated. Finally, other 
vacuum polarisation effects around flux tubes in different variants of QED, in three and 
four dimensions are mentioned. 
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1 Introduction 



Three dimensional quantum electrodynamics (QED3), with massive electrons in an irre- 
ducible two-dimensional representation, is peculiar compared to QED 4 in that it is not 
invariant under a discrete parity transformation. This breaking of parity by the massive 
fermions has a striking consequence that radiative corrections to the photon propaga- 
tor generate an additional gauge-invariant but parity-odd contribution to the effective 
lagrangian for the gauge fields M. The last is the celebrated Chern-Simons term 0. 

The two-component nature of the fermions, the breaking of parity by fermion masses, 
and the Chern-Simons term, form a triad of interlinked facts which are responsible for 
various interesting effects in QED3, one of which is the induction of fractional charge 
in the presence of a background magnetic field |], [J. To fix the notation, consider the 
lagrangian of three-dimensional QED, 

C + = ^ + (^ m7 m - ™)*+ , (1.1) 

where = <9 M + ieA^, 7 = a 3 , 7 1 = ia 1 , 7 2 = io 2 and where are two component 
spinors. (The subscript "+" in the above equations has been inserted for latter conve- 
nience.) As mentioned, for m 7^ 0, C + is not invariant Q under the 2 + 1 dimensional 
parity operation (x,y) — ► (— x,y),^ + — > a 1 ^ + . The mass parameter m in fll.lj) is al- 
lowed to take either sign (±) so as to concisely describe the two inequivalent irreducible 
representations of QED3. If the potential decribes an external magnetic field with 
dimensionless flux 

F = — I A -dr. 

then the induced charge is given by 

Q = |/d 2 *<[*+,7°*+]> (1.2) 
_£ f d 2 x £sign(£)4(x)^(z) (1-3) 

E 

= -sign(m)^ . (1.4) 

In Eq.( |1.3|) , ipE are eigenstates of the Dirac Hamiltonian H + , and so the induced charge 
is related [|l], [3| to the spectral asymmetry of H + . The last equation shows that in a 
time-independent external magnetic field, the asymmetry in the spectrum is an invariant 
depending only on the flux. There are numerous derivations of (1.4) in the literature for 



generic magnetic field configurations [I], ^ and a physical picture is as follows 0: free 
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spinors have their spin correlated with their energy, s = sl SW mg ) ; so a perturbing 
magnetic field, because of the spin-field interaction, polarises the virtual particle states 
differently according to E > or E < thus causing an induced charged cloud around 
the field. 

There are other induced quantum numbers in a time-independent external magnetic 
field, like angular momentum || J and spin || |7j S which depend only on the flux, 

F 2 

J = -sign(m) — , (1.5) 
\F\ 

S = — sign(m) , for m — > 0. (1.6) 

All of the above induced quantities depend on sign(m) which indicates the crucial de- 
pendence of the results on the physically different representations. Also, as required, 
Q is odd under charge conjugation (C) while J and S are C-even in agreement with 
the transformation properties of the corresponding single particle operators. Actually, 
in explicit calculations, the correct C-symmetry is ensured only if one computes properly 
(a)symmetrised expressions and is careful to choose regularisation schemes which preserve 
the (a)symmetry. Note also the very important point that while the expressions (1.4) and 
( |1.5| ) for Q and J are valid in the massive theory, S is an invariant (depending only on 
the flux) only in the massless limit of the theory. 

An interesting external field configuration is an Aharonov-Bohm || (AB) string for 
which the magnetic field is nonvanishing only at one point. Nevertheless, even in this 
singular case the relations (|1.4| - |1.5|) are still valid as has been verified by explicit calcula- 
tions [|l(], [11, 12, [l]|. In Ref.fllT] it was stressed that, because of the nonperiodic (in flux) 



induced charge localised around an AB string in 2 + 1 dimensions, an otherwise quantum 
mechanically unobservable string can reveal its presence to a probe electron in a Gedanken 
scattering experiment. This may be seen as an example in which quantum field (virtual 
many body) polarisation effects invalidate quantum mechanical expectations. 

While the induced charge and angular momentum have been computed for an AB 
configuration, to my knowledge the induced spin has not. Other than explicitly verifying 



the relation ( |1.6| ) for the singular AB configuration (see Sect. 2), and thereby reinforcing 
the last statement in the previous paragraph, the calculation involves a complication in 
that the limit m — > in (|1.6| ) must be taken before doing the space-integral over the spin- 
density. Actually, for an AB configuration, the only free dimensionfull parameter in the 
problem is the mass m since the magnetic field enters only through the dimensionless flux 
F. Thus for the AB problem, T = / d 2 x{^ + ^^ + ) is independent of the magnitude of m 
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and so trivially T AB = 7^ . But since T is not related in any simple way to the spectral 
asymmetry, one does not expect, in general, its massless limit to coincide with S (that 
is, the massless limit required in (|1.6Q need not commute with doing the space- integral), 
and this is discussed further in Appendix A. 

In QED^, if the electrons are massless, the theory has a chiral symmetry. This is not 
true for the lagrangian ( |1 . 1|) of QED 3 because there is no matrix which anti-commutes 
with the 7^'s. However by enlarging the representation space of the spinors in (|1 . 1| ) , one 
can also discuss chiral symmetry in massless QED 3 . Consider thus the three space-time 
dimensional lagrangian given by [|L~f 



£ = m(iD^ T^-m)^ , (1.7) 
where the /emr-component spinor \& forms a reducible representation of the Dirac algebra 

with 7 = a 3 , 7 1 = ia 1 , r j 2 = ia 2 as before. When m = 0, this lagrangian is invariant 
under a U{2) symmetry generated by J,r 5 , and — iT 3 , where F 3 = ( . j^) and T 5 = 



,i 0, 

«r°r 1 r 2 r 3 . This symmetry is sometimes referred to as a flavour symmetry. Alternatively, 
since the symmetry exists only for the massless theory, one may also call it a chiral 
symmetry. Furthermore, by writing \I/ = y ^ + ^ , with ty a (a = ±) two-component 

spinors, the Lagrangian density (|1.7|) is invariant (for any m) under the generalised 
parity operation (x,y) — > (—x,y), — > cr l ty_ and 1 i r _ — > . 

As ( |1.7| ) has symmetries similar to those of physically interesting, but technically 
difficult, four-dimensional theories like QCD, it has been used as a toy model to study 
such issues as chiral symmetry breaking. Most of the last work has been on dynamical 
chiral symmetry breaking, which involves the difficult problem of finding self-consistent 
nonperturbative solutions to the Schwinger- Dyson equations flj]]. A simpler, but still 
informative, task is to study chiral symmetry breaking induced by external electromag- 
netic fields. In Ref . fp~5f| it was found that a space-time independent external magnetic 
field yielded a nonzero value for the chiral order parameter, ( 1 l r + \E' + ) m ^ . This result 
was generalised to inhomogeneous (but still static) fields in Ref.|16| by observing that in 
an external field (that is ignoring virtual photon corrections), the chiral condensate of 
the lagrangian C is related to the induced spin density of the lagrangian C + . In terms of 
two-component spinors, Eq. ( |1 . 7| ) may be written as 

C = ^^E^^^f-am)^, (1.9) 
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so that if is an external field, ( |1.9| ) describes two decoupled systems £±. The fermion 
condensate is therefore 



C(x;m) = (y(x)*{x)) = ]T a (V a (x)V a (x)) a (1.10) 



a=± 



a C a (x;m) (1-H) 



a=± 



^aC + (x;am), (1-12) 



a=± 



where + denotes the expectation value of 1 I , + \E' + in the C + subsystem. Thus in 

an external field the properties of the condensate C(x;m — > 0) are determined completely 
by the induced spin density \ C + = in the subsystem C + which is a better 

studied problem. An immediate consequence of (|1.12|) is C(x; —m) = —C(x;m), so that 



the condensate C(x;m — > 0) can be nonzero only if the massless limit is discontinuous, 
that is it depends on m —>■ ± . That this is indeed the case was noted in Ref.|T6| and 



Sect. (3) of this paper extends the discussion of chiral- symmetry breaking given there. 



The plan of the paper is as follows. In Section 2.1 a derivation of the relation (|1.6|) 



is reviewed and in Sect. 2.2 it is verified by an explicit calculation in the AB system. 



Section 3 repeats some of the discussion on chiral symmetry breaking given in Ref. 
and supplements it with some qualitative comments. The concluding section of this 
paper summarises some of the discussion and highlights other vacuum polarisation effects 
in QED% and QED4. The appendices contain some details of the calculations presented 
in Sect. 2. Appendix B also presents a novel derivation of (|1.4|) for the AB case in the 



massless limit. The reader is informed that for continuity of discussion, some material 



from the short report jUl has been reproduced in this article. 



2 Induced Spin 

2.1 General 

Consider the induced vacuum spin in the theory defined by (|1 . 1|) , due to vacuum polari- 
sation in a time-independent external magnetic field, 

2 nwO 2 

Since m can be of either sign, the limit in ( |2.1| ) should be understood correspondingly 
as m — > 0^. On expanding the fermion operators in terms of the eigenstates ipE of the 
Dirac Hamiltonian H + = 7°7*(p i — eA l ) + my (in the A = gauge), symmetrising 
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with respect to charge-conjugation, and subtracting the infinite vacuum contribution, one 
obtains 

\c + {x-m^Q) = —lim VsignWVV's IJ". (2-2) 

The summation in ( |2.2j ) symbolically denotes a sum over discrete states and an integral 
over the continuum. In the m — > limit, the E > and E < eigenstates are related 
17| by ip-E = cr^tpE, so that only the zero modes contribute to (2.2). For F > the zero 



modes occur at E = m — > and are of the form ip ~ f q J ' wm ^ e f° r < they occur 

at = — m — > and are of the form ip ~ ( ^ ) ' ^ )eno ^ e ^ ^ ne smanes t integer greater 
than or equal to \F\ — 1; then for a given m = 0^, there are / discrete (normalisable to 
unity) states in addition to the continuum (scattering) states [T7], Eq. (|2.2j) therefore 
becomes 

1 n ( n±^ sign(m) , 

- C + (x;m -> 0*) = 2^ WbIo ■ ( 2 - 3 ) 

/ 4 £=m-*0± 



The discrete zero modes in the sum Q2.3p are localised [F7, Q around the magnetic field and 
so give a nonzero local contribution. On the other hand the continuum zero modes, being 
scattering states, give a negligible local contribution @ unless the flux is infinite. However 
the continuum states are important to establish the index theorem which is obtained by 
integrating (|2.3| ) over all space : The discrete states contribute an amount — sl S n ( m ) 
while a careful analysis shows that the scattering states contribute — S1 & n(m ) (jF\—I), 

f , 9 C + (x;m— > O 1 * 1 ) sign(m) . . . 
S = / d 2 x +V ' = \F\ . (2.4) 
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The discussion above has been for regular field configurations. If the flux contracts to 
a single point to form an AB flux string then no discrete normalisable states form even 
for \F\ > 1 because of the highly singular nature of the field. Rather, if one starts from a 
finite size tube and shrinks it to zero, the would be discrete states collapse to point-like 



states located on the string. These hidden states maintain the relations (1^, 1.5 , |1.6| ) in 
the AB case when \F\ > 1 but then they cause a breakdown in periodicity, with respect 
to the flux, of induced quantum numbers in that configuration [11], [12| . 

2.2 The AB case 

The purpose here is to verify (|2.4| ) explicitly for an AB string with \F\ < 1. The eigenstates 
of the AB hamiltonian in polar coordinates, and in the Aq = 0, e A = — 6 gauge, are 
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readily determined, 



1>. 



E,n 



(E — m)k 
4nE 



gi(n— 1)9 



E—m 



Jr 



(At) 



And 



(2.5) 



I = n — F, n = integer, \F\ < 1, F^O , 

E = ±Vk 2 + m 2 , < A; < oo, m^O, 

e n = +1 for n > 1 and e n — — 1 for n < . 

In the above, the sign of e ra has been determined by the usual requirement of square- 
normalisability |18[ of the wavefunctions and a possible ambiguity for the critical partial 



wave < l c < 1 is resolved by considering the string to be a limiting case of a flux tube of 
finite size 0. Notice that the critical partial wave spinor contains one component which 
is singular at the origin. Also note that the values F = 0, ±1, which have been excluded 
from (|2.5| ) for technical convenience, can be obtained in matrix elements (the induced 
charge, spin etc.) by taking the limits F — > 0, ±1. 

Substituting (Of) into the right-hand-side of fl2.2p and integrating over all space gives 



S = — [d 2 x fmV r dkk {Jl a „ + Jl m | |^ (2.6) 

where the massless limit has to be taken inside the space-integral. Before evaluating this 
quantity, consider first the same expression but with the massless limit taken after doing 
the spatial integral. Calling this last object 7^ , one has 

m f ,n f°° dkk 



where 

G n (kr) = J 2 n „ F (kr) + J 2 n+ Akr) + t F+ iW + • (2-8) 

The quantity T AB is discussed at length in Appendix A. 

Let us now return to the evaluation of ( |2.6j ), that is, of S. For this we need the massless 
limit of the wavefunctions in (|2.5|) . It is easy to verify that for m = 0, V-fc,n — —?Pk,n 
and if this is used indiscrimately in (|2.2| ) one obtains zero for the induced spin. The 
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resolution of this paradox lies in recalling that there ought to be unpaired continuum zero 
modes (E = ± ). In order to isolate these unambiguously one has to start with m ^ 0. 
Now when m ^ (that is before the massless limit), the states with energy E ~ ±m 
correspond to values of k satisfying k <C \m\. Thus one considers the A — > limit in 
( |2.5| ) to identify the nonvanishing threshold states. In this limit only the critical partial 
wave survives: For 1 > F > its limiting form is 



while for > F > — 1 it is 



k_ fJ- F {kr) 
2tt V 



lim \l/„ c =i = sign(m) \l — ( ' \ , E fa m, F > . 



(2.9) 



lim * 

fc-^O 



n c =0 




£? w -m, F < . 



(2.10) 



(The Bessel function Ji_pi(A;r) has not been expanded in k because kr can be large even if 
k is small). When finally the limit m — > is taken, these ( |2.9| - |2.1U| ) will be the only states 
which will contribute to the spin, the contribution from the \E\ > states cancelling in 



pairs as mentioned above. Thus for \F\ < 1 the [19] scattering zero modes (2. £- 2. 10 ) 
contribute to (|2.2|) the amount 



sign(m) r x dk 



o 2tt 



k Jl ]F \(kr)\* 



(2.11) 



A 2 sign(m) 
IQir 



{j 2 |F| (Ar)- J-\ F \-i(Xr) J HF , +1 (Ar)}| F , (2.12) 



where A |m|. Note that taking the limit |m| — > (and so A — > 0) on the right-hand- 
side of ( |2.12j ) results in a null contribution to the (local) spin-density in accordance with 
the remarks in the previous subsection that for \F\ < 1 only the dilute continuum zero 
modes contribute to the density. However the space-integral of the density in ( |2.12| ) is 
non-negligible and is given by 



27rsign(m) / — r(^ + ^+) m ^ 
Jo 2 



XH(0) F 

lo ' 



A 2 



J dr r | J_|ir|_i(Ar) J_ 



1 r°° 
8 Jo 

i(-2|F|) 

\F\ 



x { J_|f|_i(a;) J_| F |+i(x)} [^ 



(^)}| F 
(2.13) 
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and so 

2tt / dr -0Mr + ) m ^ o = -sign(m)i-l, |F| < 1 , (2.14) 
jo 2 4 

which agress with the index theorem (|1.6|) . Note that the integral in (|2.13| ) is independent 
of the parameter A so that the final result is not sensitive to the slight ambiguity in the 
choice of this cutoff. The evaluation of the integral (|2.13| ) is described in Appendix B. 

To summarise the results of this section, the induced spin in the AB system was 
computed for \F\ < 1 and the index formula ( |1.6| ) verified. It was shown that isolating 
the contribution of the threshold states required a careful consideration of the massless 
limit. For \F\ < 1 the contribution to S is smeared over all space due to the nonlocalised 
nature of the continuum (scattering) threshold states and thus the spin-density (in the 
required massless limit) is zero. For \F\ > 1 there are also discrete states for any finite size 
flux tube which give a localised nonzero contribution to S (see the analogous discussion 



for the charge in (TT|, 0). As mentioned earlier, these discrete states collapse to point-like 
states in the limit of a zero-radius flux tube. 

The presence of induced quantum numbers around the QED 3 AB string implies that 
it is detectable (through a scattering experiment) even for integral flux jrij. However, 
there is a return flux at infinity (actually, at r ~ where /i ~ e 2 is the dynamically 
generated photon mass) complementary to the AB flux in the string, and one can ask if 
the vacuum polarisation effects associated with that return flux can aid in restoring the 
invisibility of the integral flux string. If the mass m of the electron is much larger than 
/i, then for any F the polarisation effects (|1.4| - |1.5| ) are localised within 1/m of the string 
and are probably not significantly JlT] affected by the return flux effects located at 1/fi 
. For F exactly an integer the induced effects are exactly on the string and so should be 
even more immune from the return flux effects. 

The situation is more complicated for massless electrons, which is the limit in which 
the induced spin (|1.6|) is defined. However even in the massless case, the contribution 
from the collapsed threshold states (for \F\ > 1 ) is highly localised so it is not very clear 



to me if these bound states will be completely nullified in their contribution to ( |1.4[ - |1.6| ) 
by feedback effects. One does however expect the contribution from the scattering states 
to be cancelled by feedback . 

There are other subtleties associated with the massless theory. For example, the factor 
of sign(m) in ( |2.14|) , if averaged over, suggests no polarisation effects in the massless limit. 
However if one adopts the point of view that the theory is sensibly defined only for massive 
electrons (to define an unambiguos zero of energy) then no matter how small that mass is 
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with respect to e 2 , the sign(m) factor distinguishes between the two inequivalent physical 
representations one is working in and so should not be averaged over. Indeed, even if one 
starts with m = in (]1.7|) , a small nonperturbative mass is generated dynamically [14 ]. 



This suggests that the ambiguities of a truly massless theory are probably academic and 
one should at most speak of a massless limit, where the fermion mass is very much smaller 
than the other dimensionfull parameter, e 2 . 

The discussion in this section has been for the parity non-invariant lagrangian ( |1 . 1| ) . 
In the parity-invariant theory ( |1.7| - 11.9| ) there is no photon mass generation due to cancel- 
lations between the contributions of the "+" and "-" spinors, and so no induced charge 
like (ITT 



3 Chiral symmetry breaking 

Let q be the Dirac fermion operator in four space-time dimensions. Then to study dy- 
namical chiral symmetry breaking, one looks at the condensate < qq > in the massless 
limit of the full dynamical theory. In general this is difficult to do analytically and an 
intermediate step that can give some insight is to look at the condensate for a particular 
background gauge-field. In Ref.[^l|] it was shown that if < qq > A is the space-time average 
of the condensate in the background field A^, and p.o(A) is the density of eigenvalues of 
the background operator ilp per unit space-time volume, then 

<qq> A = np D (0). (3.1) 

That is, chiral symmetry breaking in a given background field is related to the spectral 
density at zero energy of the operator i p. If one can parametrize the right-hand-side 
of ( |3.1| ) explicitly in terms of the background field A^, then by integrating (in the path- 
integral) Eq.( |3.1| ) over all possible gauge-field configurations, one can in principle study 
dynamical chiral symmetry breaking in the full theory. In practice however, it is already 
quite difficult to achieve the first step above : that is, to make statements about ( |3.1| ) for 
arbitrary configurations. 

However for three-dimensional QED described by the parity-invariant lagrangian ( |1.7|) , 
much can be said. For external time-independent magnetic fields with flux F, the con- 
densate is given by (|1.12| - |2T3| ) to be 

Cforo-O*) = -sign(m) ]T ^ E \v , (3-2) 

E=m-^0 ± 

where ipE are the eigenstates discussed in the previous section. It is quite clear that (|3.2|) 
is just a special case of ( |3.1| ), obtained by restricting to time- independent magnetic fields. 
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However ( |3.2| ) contains some unique features which enhance its usefulness in spite of the 
restrictions. For any finite flux, the spectral density on the right-hand-side of ( |3.2|) is 
dominated locally by a finite number of discrete (normalisable to unity) zero modes. The 
number of such modes is determined solely by the flux as discussed earlier. Furthermore, 
the integrated density also depends only on the flux, 

POO _ 

/ <Px(* + * + ) m ^ Q = - sign(m) |F| . (3.3) 
Jo 

This last is significant for it means that for any configuration with nonzero net flux, the 
chiral condensate must be nonzero at some place. Thus one can decide if chiral-symmetry 
is broken or not by just knowing the flux rather than the details of the magnetic field. 
This also suggests that the contribution of static magnetic fields to the full dynami- 
cal condensate is not negligible (of nonzero measure in the path-integral) and makes it 
plausible that dynamical chiral symmetry happens in the full dynamical theory (that is, 
without external fields). Of course the last is already known to happen (JTJ[] and so the 
real challenge is to see if one can use ( |3.2| - |3.3| ), integrated over arbitary static magnetic 
fields (of variable fluxes) to get an estimate for the dynamical condensate. By ignoring 
time-dependent fields, and in particular electric fields, there is no a priori reason such an 
estimate would be close to the results of Ref.]14|], but still the exercise might be revealing 
if doable. (This will not be attempted here.) 

Another important feature of ( |3.2| ) and ( |3.3| ) is the sign(m) factor. As noted in Ref . [|T6[ 
this means that the direction of chiral symmetry breaking in external magnetic fields is 
correlated with the infinitesimal breaking of that symmetry brought about by a mass 
m — (see the discussion of massless theories near the end of the last section). 

Let me mention now some perturbations which could destroy the condensate ( |3.2| - |3.3| ) 
formed by static magnetic fields. Viewing the condensate as an ordering of spins in a 
magnetic field, one expects a superimposed electric field to reduce such an ordering and 
so cause a reduction in the magnitude of the condensate. At least if the electric field is 
weak then one might be able to check this in a perturbative treatment as was done in 



Ref.p2| for a Chern-Simons theory. 

Thermal fluctuations are also naturally expected to disrupt any ordering (symmetry 
breaking). For example, the charge ( |1.4|) has been computed @, |23| for slowly varying 
magnetic fields at nonzero temperature and density, showing a rapid "washing out" at 
high temperature. More remarkably, for any non-zero temperature, the charge was shown 
to vanish in the massless limit; a simple explanation of this given in |3j was that in the 
massless theory the charge is due to the zero-energy modes which are susceptible to large 
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fluctuations in their occupancy at nonzero temperature. Since the chiral condensates 
( |3.2f - |3.3D also receive support only from the zero modes, the qualitative argument given 



above suggests that they should also vanish at nonzero temperature. This has recently 



been noted in an explicit calculation for a constant magnetic field p4 |. 

4 Conclusion 

In this paper were studied some phenomena in three-dimensional QED in an external 
magnetic field. The breaking of chiral symmetry by an external magnetic field in the 
reducible represent at iuon of QED 3 (|1.7|) is equivalent to the polarisation of spin in the 
irreducible theory ( |1.1| ). Thus chiral symmetry breaking in the toy-model ( |1.7|) may be 



heuristically understood as an ordering of spins in a magnetic field. Then any effect which 
interferes with this ordering, such as an electric field or thermal fluctuations, would be 
expected to diminish the breaking of chiral symmetry. 

Quantitatively, the chiral condensate in an external magnetic field is controlled by 
the total flux of that field ( |3.2| ) as it determines the number and nature (discrete and 
localised versus continuum and non-localised) of the zero modes of the massless Dirac 
equation. For any finite flux, the localised zero modes give the essential non-negligible 
contribution to the chiral-condensate. The space-integral of the condensate receives a 
nontrivial contribution also from the dilute scattering zero modes and the net result is 



an invariant depending only on the flux ( |3.3|) . This index formula gives the phenomena a 



certian robustness and might be useful in obtaining estimates of the chiral condensate in 
the full dynamical theory as described in Sect. 3. 

The polarisation effects (|1.4j , |1.5| , |1.6| ) in the theory (b_l) are also invariants depending 



only on the flux. It has been verified ( |]IU| , O, |T^, [T3| and Sect. 2) that the formulae ( |1.4| , 
|1.5| , |1.6| ) hold also for singular configurations such as the AB flux string. Consequently, 
as a result of the localised threshold states which form every time the flux passses an 
integer, there accumulate around an AB string vacuum polarisation effects which are not 
periodic in the flux. This exposes an integral flux (F ^ 0) string to external probes 
[TT| . By contrast, such polarisation effects cancel for flux tubes in QED4 due to the four 



component nature of the fermions which preserve parity TQ. However this does not 



rule out other polarisation effects (as yet undiscovered) which might also make an integral 
AB flux string in QED4 visible. 

In the model ( |1.7| ) the net induced charge around a magnetic field vanishes because 
the contribution from the " + " and " — " two-component spinors cancel (since they 
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have opposite signs of mass). One can consider a modification of ( |L7| ) which is also 
parity invariant but in which the net induced charge is nonzero. This model is obtained 
by replacing the charge e in ( |1.7| ) by the matrix e* = e ^/V This means that 

the ± spinors now have opposite charges in addition to oppositely signed masses and 
their respective contribution to the induced charge will be additive. This last model 
has actually been used as a parity-invariant effective lagrangian for high-temperature 
superconductivity, but the gauge-field which couples to the charge e* there is a "statistical 
gauge field" rather than the true electromagnetic field. 

There certainly exist vacuum polarisation effects around flux tubes, even in parity-odd 
three-dimensional QED ( |f ,f|) , which are not caused by parity breaking effects like fermion 
masses. For example, there are vacuum currents which circulate around thin magnetic flux 
tubes and are due to the AB effect acting on virtual electron-positron pairs |26|, [11], |12|| . 
These currents occur also for scalar theories and in 3 + 1 dimensions. A possible role 
for them has been proposed recently in a picture of a nonperturbative vacuum of QED 4 



27, 28 



Whether one uses three-dimensional QED as a simplified model in particle physics, as 
an effective theory in condensed matter physics, or as a field theory of intrinsic interest, 
it is clear that much more remains to be understood. The focus in this paper has been on 
extending previous discussions of polarisation effects due to magnetic fields. The inter- 
ested reader is referred to some other recent works [2|| in which magnetic fields are also 
employed to unravel the dynamics of QED3. 
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Appendix A 

Here the quantity T AB mentioned in Sect. 2 is evaluated, and some subtleties involved 
are discussed by comparing the calculation with the corresponding one for the induced 
charge Q. 

Simplifying for the AB case one obtains, 

r°° dk k 



-1 

47T 



m 



E 



rn 



o y/k 2 + m 2 
00 dkk 



{ J e n (l-l) + J L(l)} \o 



4tt Jo y/k 2 + 



n>0 



where 



G n (kr) = J 2 n „ F {kr) + Jl +F {kr) + Jl„ F+1 {kr) + Jl +F+1 {ki 



Note the manifest C-symmetry, F — > —F, of the expression for G n . By writing 

1 

y/k 2 + m 2 \fn Jo 
the k— integral in ( |A.2|) may be performed using [30 



4= f°°da e ~^ k2+m2 \ 



/ e~ a k Jl(kr) kdk = — 
Jo mV 7 2a 



2^ e ~^Mi 



to give 



(*+*+>- 



m 



da 



Aix^fn Jo a 2 



-a m —z 



»n>0 



where 



H n {z) = I n -.p(z) + I n+ p{z) + I n -p + i{z) + I n+ p + \(z) 



(A.l) 
(A.2) 

(A.3) 
(A.4) 

(A.5) 



(A.6) 



(A.7) 



and z = r 2 /2a 2 . Now consider the spatial / rdrdO integral over the density in ( |A.6| ). The 
a— integral can then be done explicitly because it reduces to the simple Gaussian 



da e 



7T 



2 m 



(A.8) 



after the change of variable r 2 = 2a 2 t in the r— integral. Thus one is left with 



d 2 x(^ + y + )- 



sign(m 



(A.9) 
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Assuming that the order of integration in ( |A.9| ) may be interchanged with the summation, 
the t— integral may be evaluated by taking appropriate limits of the formula 



L 



oo KP+ 
at- " 



e- at I p+u (bt) = - . (A. 10) 



For example, for any integer n, 

i 

e-\l n+F (t)-I n (t)) = -F. (A.ll) 
Thus the final result for the space-integral of (|A.1|) is 

'd 2 x(^ + ^ + ) + = - Sign ] (m) J2 (F-F + F-F) = 0. (A.12) 



n>0 



Therefore T AB = |(Eq.([AT2])) = 0, if the space- integral above may be performed before 
the infinite partial wave summation. 

By way of comparison consider the charge (|1.3|) evaluated for the AB case with | F | < 1 . 
Substituing ( |2.5| ) into ( |1.3| ) and simplifying gives [PH] , 



—me sign(F) f f f°° dk k 
eF 



4, j hi: whi^-^) 



= — sign(m)^-. (A. 15) 

Apart from a factor 2e, the expression (|A.13|) differs from T AB by a relative minus sign 
between the two Bessel functions in curly brackets. Consequently the sum over partial 
waves in ( |A.13|) results in several cancellations to give ( |A.14j ). By writing J\p\ — J-\f\ — 
(J\ F \ - Jo) - (J-\f\ - Jo), Eq. (|A.14|) may be evaluated using (|AT^ , |AT5| , |A.11|) to get (|A.15|) . 
(This computation is slightly different from the one in Ref. JTT| where ( |A.11| ) was not 
used) . 

Since Q is C-odd no F = subtraction is necessary in (|A.13|) as that piece is auto- 
matically zero. On the other hand since T (and S) are C-even, a F = subtraction is 
necessary to remove an infinite bare vacuum contribution. The last infinity is due to the 
fact that the spin of free fermions is given by s = \ sign (mE) so that the spin of the bare 
ground state, J2e sign(E) s ~ sign(m), is divergent. 
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If the inifinite summation over partial waves in ( A.13 ) is regulated by a cut-off. 



n 



max = N, (N > 0) which is later taken to infinity, one obtains (see |20| 



_ -me f rf U E f " (A,,-^,,) (A.16) 

= ^ & / a /; -j^= ({ ji f - 4} + {j>_ - 4 +1+f }) (A.i7) 

It is easy to verify using the formulae given earlier that if the limit N — > oo in (|A.17|) is 
taken after the space- integral, the answer is zero as the contribution from the iV-dependent 
Bessel functions cancels the contribution from the first two. However, as noted in Ref. ||20|[ . 



in the large N limit the contribution from the N- dependent pieces recedes to infinity and 
it may be interpreted as representing the oppositely charged polarisation cloud "repelled" 
by the localised flux. So if one is interested in the local polarisation effects, the N — > oo 
limit should be taken before the space- integral in (|A.17|) . Note also that for any finite 
iV the expression ( |A.17 ) is no longer C-odd, this asymmetry being broken by the N- 



dependent pieces. This may be remedied by taking the C-odd part of the truncated sum 
in (|A.16|) . The rest of the discussion in this paragraph still goes through then. 

Returning now to the evaluation of T AB , consider regulating the partial- wave sum 
(A..9) in a C-even way as for the charge above. Clearly, the result is still zero if the 



space-integral is done first before the summation. This suggests, by analogy with the 
charge case above, that in this order (space-integral before summation) of calculating T 
one is adding the local polarisation effects to the opposite non-local ones. Unfortunately, 
because there are no pair-wise cancellations among the partial waves in (|A.9|) as in the 
charge case, it is not obvious which of the partial waves are contributing to the local po- 
larisation effects and which pieces may be interpreted as the opposite effects which recede 
to infinity in the N — > oo limit. Fortunately however this question, though interesting, 
is not of direct relevance because the physical object S as defined in ( |2.4j ) differs from T 
in that the massless limit has to be taken before the space-integral (as suggested by the 
general derivation sketched in Sect. 2 ). The evaluation of S is described in the main text. 
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Appendix B 

Here the integral in Eq.( |2.13| ) is evaluated. It is 

/ dx X \ J-\F\-l{x) J-\f\+i(x) — J-\{x) Ji{x)\ 

J 

dx X \ J-\F\-l(x) J-\f\+i(x) + Ji(x)j 



(B.l) 



with |F| < 1. The integral ( |B.1|) is well defined but inconvenient to evaluate as a 
whole. Therefore the integral over the F-dependent and F-independent Bessel functions 
in (B.l| ) will be performed separately. This will require an intermediate regularisation. 
The method adopted here is the following. One has an integral |3(| 



Mt) Ju{t) r a dt 



with 



r(a) T( v+li - a+1 ] 



u + fi + l>a>0 



(B.2) 



(B.3) 



By setting fi = —\F\ — 1 and v = —\F\ + 1 in ( B.2 ) and choosing F and a so that 
the condition ( B.3| ) is satisfied, the F dependent integral in ( B.l|) may be performed by 
analytically continuing the right-hand-side of (|B.2|) to a = — 1 and the full range of F. 



The F-independent integral in (|B.1|) may be evaluated similarly. The net result for ( B.l|) 
is -2|F|. 

It is useful to have a cross-check on the above manipulations with ( |B.2| ). For this 
purpose, consider the following calculation of the charge fll.3| ), for the AB system with 
flux \F\ < 1, in the massless limit of the theory. Using the continuum threshold states 
given by (Ol -PTTDf) and the mapping ip E — a^ipE between the non-threshold states, one 



obtains (cf. the calculation of the spin in 2.12 - 2.14 ) 

lim Q = — sign(m) sign(F) lim / 2irrdr / — J-\F\(kr) \q (B.4) 

2 A->0 Jo JO 27T 



— sign(m) sign(F) \F\ 
2 



(B.5) 



-sign(m) F 



:b.6) 



where A <C \m\. In (|B.4|) a F = subtraction has been inserted to remove divergent 
F-independent pieces. Such a subtraction was not required in the calculation of Q (for 
m 0) by the conventional method in Appendix A (see |A.13| ) and so presumably it is 
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required here because of the ambiguity in the cut-off A. The integral in ( |B.4| ) has been 
evaluated using (p.2|) and the final result ( |B.6D agrees with the expected formula Q1.4J) 
which holds independent of the magnitude of m. (Please note however that the charge- 
density is not a topological invariant and it also depends nontrivially on the mass m. 
Thus, for example, for \F\ < 1 the charge-density is infinitely dilute in the m — > limit 
(receiving contributions only from continuum zero modes (|B.4j )) but is concentrated within 
a distance 1/m of the localised field for m 7^ (in which case it receives contributions 
also from non-threshold continuum states) p], |12|.) 
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